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The standard quantum search lacks a feature, enjoyed by many classical algorithms, of 
having a fixed point, i.e. monotonic convergence towards the solution. Recently a fixed 
point quantum search algorithm has been discovered, referred to as the Phase- it/ 3 search 
algorithm, which gets around this limitation. While searching a database for a target 
state, this algorithm reduces the error probability from e to e 2q+1 using q oracle queries, 
which has since been proved to be asymptotically optimal. A different algorithm is pre- 
sented here, which has the same worst-case behavior as the Phase- n/3 search algorithm 
but much better average-case behavior. Furthermore the new algorithm gives e 2q ^~ 1 
convergence for all integral q, whereas the Phase- n/3 search algorithm requires q to be 
(3 n — l)/2 with n a positive integer. In the new algorithm, the operations are controlled 
by two ancilla qubits, and fixed point behavior is achieved by irreversible measurement 
operations applied to these ancillas. It is an example of how measurement can allow us 
to bypass some restrictions imposed by unitarity on quantum computing. 
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1 Introduction 

Quantum computing gives us a powerful computational framework, by exploiting the super- 
position and entanglement phenomena exhibited by quantum systems. A famous example 
of this power is the quantum search algorithm , which provides a quadratic speedup over 
classical search algorithms. This quantum search algorithm consists of an iterative sequence 
of selective inversion and diffusion type operations. Each iteration results in a fixed rotation 
(which is a function of the initial error probability) of the quantum state in a two-dimensional 
Hilbcrt space formed by the source and the target states. If we choose the right number 
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of iteration steps, we stop very close to the target state, else we keep on going round and 
round in the two-dimensional Hilbert space. To perform optimally, therefore, we need to 
know the right number of iteration steps, which depends upon the initial error probability or 
equivalently the fraction of target states in the database. We can estimate this fraction to 
the desired accuracy using various "amplitude estimation" algorithms but that requires 

additional queries. When the total number of queries is large, the additional queries do not 
matter much, but when the total number of queries is small, the overhead can be unacceptably 
large. 

In this paper, we address the problem of finding an optimal quantum search algorithm in 
situations where, (i) we do not know the initial error probability (perhaps only its distribution 
or a bound is known) , and (ii) the expected number of queries is small (so that every additional 
query is a substantial overhead). Such situations occur in pattern recognition and image 
analysis problems (where each query is a significant cost), and in problems of error correction 
and associative memory recall (where the initial error probability is small but unknown). In 
these cases, the quantum search algorithm and its "amplitude estimation" enhancements are 
not of much use. We look for different variations of the quantum search algorithm that ensure 
amplitude enhancement, and also outperform classical search algorithms. 

The strategy, familiar from classical computation, is to construct an algorithm that "con- 
verges" towards the target state. That is impossible to do by iterating a non-trivial unitary 
transformation; eigenvalues are of magnitude 1, so the best that can be achieved under such 
conditions is a "limit cycle" and not a "fixed point" . As described above, this is indeed what 
happens in case of the quantum search algorithm. To obtain an algorithm that converges 
towards a fixed point, some new ingredient is required, and several possibilities come to mind: 

(a) Some property of the current computational state offers an estimate of the distance to 
the target state. This estimate can be used as a parameter to control the extent of the next 
iterative transformation, e.g. the Newton-Raphson method to find the zeroes of a function. If 
such an estimate is not available, or is too expensive, we have to look for a different method. 

(b) Suitably designed but distinct operations are performed at successive iterations. Such a 
method can converge towards a fixed point even using unitary transformations, as exemplified 
by the recently proposed Phase- ir/3 search algorithm 11. We will describe this in more detail 
below. 

(c) Irreversible damping is introduced in the algorithm without explicit use of any property 
of the target state. With the right type of irreversibility (i.e. when all eigenvalues of the 
fixed iterative transformation are less than 1 in magnitude), the algorithm converges, e.g. the 
Gauss-Seidel method for solving a set of linear algebraic equations. Within the framework of 
quantum computation, such an irreversibility can be introduced by projective measurement 
operations, and the algorithm of this paper falls in this category. 

By construction, fixed point algorithms provide two attractive features, which are not 
commonplace in generic quantum algorithms: 

(1) The initial state is guaranteed to evolve towards the target state, even when the algorithm 
is not run to its full completion. 

(2) Any errors due to imperfect transformations in earlier iterations are wiped out by the 
subsequent iterations, as long as the state remains in the problem defining space. 

These are motivation enough to explore fixed point quantum algorithms, even if they have 
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other limitations. 

Let us consider an unsorted database in which a fraction / of items are marked, but we 
don't have precise knowledge of /. We run a particular algorithm which has to return a single 
item from the database. If the returned item is a marked one, the algorithm has succeeded, 
otherwise it is in error. Without applying any algorithm, if we pick an item at random, then 
the probability of error is e = 1 — /. The goal of the algorithm is to minimize the error 
probability, using the smallest number of oracle queries. If / is sufficiently small, then we can 
use the optimal quantum search algorithm to obtain a marked item using 0(1/ ^/f) queries. A 
few more queries to estimate / or to fine-tune the algorithm is not a problem, because overall 
we gain a quadratic speed-up compared to the classical case requiring exhaustive search. But 
when / is large, the number of oracle queries is small, and the quantum search algorithm 
doesn't provide much advantage — in particular the rotation can overshoot the target state. 
In such a situation, a simple classical algorithm (select a random item and use a query to 
check if it is a marked one) would outperform the quantum search algorithm. The same 
considerations apply to the more general amplitude amplification algorithms |2J [3] . There 
the initial quantum state is a unitary operator U applied to a given source state \s), and the 
probability of getting a target state after measuring this initial state, |J7t s | 2 , is analogous to 
/. The probability of error, which has to be minimized using the smallest number of queries, 
is the probability of getting a non-target state after measurement, e = 1 — |J7t s | 2 - 

One of us has recently proposed a fixed point quantum search algorithm 0], which we 
refer to as the "Phase- 7r/3 search". It shows that by replacing the selective phase inversions 
in quantum search algorithm by selective 7r/3-phase shifts, the quantum state monotonically 
moves closer to the target state. Remarkably, this convergence is achieved using reversible 
unitary transformations and without ever estimating the distance of the current state from 
the target state. Explicitly, if the initial error probability after applying the operator U to the 
source state |s) is e, then the error probability after applying the operator URs^ 3 U^R^ 3 U 
to |s) becomes e 3 (here R^ 3 and R^ 3 are 7r/3-phase shift operators for the source and 
the target state respectively). Recursive application of this transformation n times, which 
requires q\ = 3qi_i + 1 (q% = 1) oracle queries at the i th level, makes the error probability 
e 3 ™. Thus e = is the fixed point of the algorithm, and the error probability decreases as 
e 2?+1 as a function of the number of queries q. This e 2q+1 performance has been shown to 
be asymptotically optimal a 5 . Note that the best classical algorithm can only decrease the 
error probability as e q+1 (not e 9 , since the last iteration does not require a query), and thus 
the Phase-7r/3 search improves the convergence rate by a factor of 2. Alternative approaches 
to the Phase-7r/3 search have been formulated 5 6, which achieve the same optimal behavior 
using different recursive schemes. An important limitation of the Phase- 7r/3 search, as well 
as its alternatives, is that the e 2q+1 performance is obtained only for a restricted number of 
oracle queries, q = (3™ — l)/2 (for positive integer n). 

Here we present a new implementation of fixed point quantum search, that gives us the 
optimal e 2q+1 performance for all positive integer values of q. It uses a new kind of quantum 
search, where the oracle and diffusion operations are controlled in a special way by two ancilla 

a Application of a general operator UR®U^ RfU to \s) changes its component along \t±) by the scale factor 
| exp(^ (8 — <f>)) — 4 sin ^ sin ^ \Uts 2 1 ■ In the asymptotic limit, \Uts\ —> 1, and the scale factor is minimized by 
the choice 8 = <f> = ^ . 
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qubits and their measurement. The same transformation is repeated at every iteration, but 
since the transformation is made non-unitary by measurement, the quantum state is able 
to monotonically converge towards the target state. Thus the algorithm is a novel example 
of how measurement can allow us to bypass restrictions imposed by unitarity in quantum 
computing. 

In the next section, we begin with a simple scheme for fixed point quantum search, and 
then modify it to present the actual algorithm. We analyze the algorithm in section 3, and 
discuss its features in section 4. 

2 Algorithm 

To obtain a fixed point quantum search, we have to find an algorithm that successively 
decreases the probability of finding a non-target state. Consider a quantum register whose 
states encode the items in the database. Let us say that the initial state of the register 
is U\s) = sin#|f) + cos#|i^) — an unspecified superposition of target \t) and non-target \t±) 
states. (Here, without loss of generality, we have absorbed arbitrary phases in the definition of 
\t) and \t±}.) To this register, we attach an ancilla bit in the initial state |0). Then we perform 
an oracle query, and flip the ancilla when the register is in a target state. Now if we measure 
the ancilla, outcome 1 tells us that we are done, and measurement of the register will give us a 
target state. Outcome tells us that the register is in the superposition \t±) of the non-target 
states. The probability of outcome is equal to the initial error probability e = cos 2 9. To 
decrease this probability, we apply the "diffusion operator" UI S U^ to the register, conditioned 
on the measurement outcome being 0. That reflects |tj_)|0) about J7|s}|0) to give the state 
sin2#|t)|0) + cos2#|i^)|0). The error probability has thus decreased by the factor cos 2 29. 
Iterating the sequence of oracle query and diffusion operations n times, the error probability 
is reduced to cos 2 6 cos 2 ™ 29. For n = 1, the error probability is e(2e — l) 2 = 4e 3 - 4e 2 + e. 
A comparison of this expression with the corresponding result e 3 of the Phase- tt/3 search is 
shown in Fig.l. We see that for e > 1/3, this simple scheme is better than the Phase- tt/3 
search, but it becomes worse if e < 1/3. 

Our goal is to find an algorithm which gives optimal convergence for all values of e, without 
knowing any bounds that e may obey. Figure 1 shows that, with one iteration of the simple 
scheme described above, the error probability monotonically drops from 1 to as e decreases 
from 1 to 1/2. That gives an intuitive idea for getting a better directed quantum search — 
somehow set a lower bound of 1/2 for e, or equivalently an upper bound of 1/2 for the fraction 
of target states /. The natural upper bound for / is 1, but we can easily make it 1/2 by using 
an extra ancilla in the state |+) = (|0) + |l))/\/2 and performing an oracle query controlled 
by this ancilla. This logic suggests the following g-iteration algorithm: 

(i) Attach two ancilla qubits in the |0) state to the source state register \s), i.e. \s) — > 
|0)|s)|0). (In what follows, we refer to the former ancilla as ancilla- 1 and the latter 
ancilla as ancilla- 2.) 

(ii) Apply H <g> U ® I to this extended register to prepare the initial state of the algorithm, 

l+>(tf|*»|o>. 

(iii) Iterate the following two steps, q times: 

Step 1: If ancilla- 1 is in the state |1), then perform an oracle query that flips ancilla-2 
when the register is in target state. 
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Figure 1 : A comparison of error probabilities after one oracle query for fixed point quantum 
search algorithms, (a) is the result for the simple scheme described at the beginning of section 
2, and (b) is the result for the Phase- 7r/3 search algorithm. 



Step 2: Measure ancilla-2. If the outcome is 1, the register is certainly in the target 
state, so exit the iteration loop. If the outcome is 0, then apply the joint diffusion 
operator (H <g> U)Io s {H ® U)t to the joint state of ancilla-1 and the register, 
(iv) After exiting or completing the iteration loop, stop the quantum algorithm and measure 
the register. 

The quantum circuit for this algorithm is shown in Fig. 2. 
3 Analysis 

Let us analyze the algorithm step by step. The initial state is 



(if 0*7® J)|0)|s)|0) = ( |0 ^^ ) (sin0|t)+cos0|t±))|O> . (1) 



The initial error probability is e = cos 2 9, and the initial success probability is / = 1 — e = 
sin 2 9. We work in the joint search space of ancilla-1 and the register, denoted by the subscript 
j, where all the states |0)|£), |0)|fj_), |l)|*x) act as non-target states. Only the state \tj) = 
acts as the target state in this joint search space, and let us represent the superposition of all 
non-target states by \tj). In the joint search space, the initial state is 

l^>=(^)l^>|0> + ^>|0>, (2) 
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Iterate q times 



ancilla-1 |0) — H 



register \s) U 



ancilla-2 |0) 




1: exit 



Figure 2: Quantum circuit for our q- iteration algorithm. 



where the unit vector \t'j) is 



T / sin cos#.„, . , cos# lH 
|t$) = N ( — |0)|t) + — |0)|t ± ) + — 11>|*±> 



V2 



V2 



V2 



(3) 



and the normalization factor N is (cos 2 9 + | sin 2 9) = \/2/(l + e). For later reference, 
note that the error probability after measuring the joint state \t'-) (i.e. the probability of 
finding the register in the non-target state |tj_)) is 



|(0|(U|i')| 2 + |(l|(t±|i')| 2 = N 2 cos 2 e = N 2 e . 



(4) 



Step 1 of the algorithm, using an oracle query, flips ancilla-2 when the joint register state is 
\tj). In step 2, we measure ancilla-2. If the outcome is 1 then we stop the algorithm, because 
the register is in the target state. The probability of getting 1 is sin 2 9/2 = f/2 — we have 
effectively put an upper bound of 1/2 on the success probability using ancilla-1. If the outcome 
is 0, which has probability 1/N 2 , then the joint state is \tj). In this case, we apply the joint 
diffusion operation using the joint source state \sj) = |0)|s). The joint diffusion operation is 
a reflection about (if® t/)|0)|s) = Uj\sj), in the two-dimensional Hilbert space orthogonally 
spanned by \tj) and \t'-) as shown in Fig. 3. The state Uj\sj) makes an angle 9j, defined by 
sin 2 9j = sin 2 9/2, with the state \t'j). So reflecting \t'j) about Uj\sj) gives us a state that makes 
an angle 29j with \t'j). Its component in ^-direction is cos 28j = 1 — 2sin 2 0j = 1 — sin 2 9 = e. 
Thus the joint diffusion operation produces the final state 



|V/> = UjLtfty = Vl^T 2 \t 3 ) + e\t' 3 ) 



(5) 



After measuring \ipf), the probability of getting |^) is e 2 , so the total probability of getting \t'j) 
after one iteration is N~ 2 e 2 . Iterating the algorithm will keep on decreasing this probability 
by a factor of e 2 at each iteration, and after q iterations it will become N~ 2 e 2q . So the net 
error probability after q iterations is (cf. Eq.(4)) 



e q = N 2 e(N- 2 e 2c >) - e 2 « +1 , 



(6) 
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Figure 3: Step- by-step quantum state evolution in our q- iteration algorithm. The double 
arrows relate the equivalent states in the original search space (left) and the joint search 
space (right). 
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which agrees exactly with the corresponding result for the Phase- 7r/3 search. 
The above analysis allows us to also deduce the following features: 

(1) e = 1 can be made a fixed point of the algorithm, instead of e = 0, by effectively 
interchanging the roles of \t) and \t±). This is achieved by flipping ancilla-2, only when 
the joint state is Then the probability of finding the register in the state \t), after q 
iterations, becomes (1 — e) 2q+1 . Note that the same behavior can be obtained in the Phase- 
7r/3 search, by replacing either R^ 3 with R t 7r ' 3 or Rs with R s ?r ^ 3 . This e = 1 fixed point 
can be useful in situations where certain target states are to be avoided, e.g. in collision 
problems. 

(2) When we use fixed point quantum search to locate the target state in a database, the 
initial error probability is e = 1 — /. The number of oracle queries required to reduce this 
probability to o(l) obeys 

(1 - ffi +l = o(l) S ^ U e -^ +1 >^ (l). (7) 

Thus we need q = 0(1/ f) oracle queries to find the target state reliably. The same is true 
of the Phase- ir/3 search as well as its alternatives This scaling of fixed point quantum 

search is clearly inferior to the 0(\/\/J) scaling of the quantum search algorithm. Still, as 
discussed in the introduction, fixed point quantum search can be useful in situations where / 
is unknown and large. 

(3) It is possible to stick to unitary operations throughout the algorithm and postpone mea- 
surement till the very end. In such a scenario, the unmeasured ancilla-2 has to control the 
diffusion operation and all the subsequent iterations (i.e. they are executed only when ancilla- 
2 is in the |0) state), and it cannot be reused in the iteration loop. We need a separate ancilla-2 
for every oracle query, to ensure that once the states \tj) and \t'j) are separated by an oracle 
query in an iteration, they are not superposed again by subsequent iterations. The whole set 
of q ancilla-2 can be measured after the iteration loop, in a sequence corresponding to the 
iteration number, to determine the target state. In this version, the unitary transformation is 
different for each iteration (because each iteration involves a different ancilla-2 and different 
controls), and the unitary iterations converge to a fixed point. 

4 Comparison 

Now we can point out some advantageous features of our algorithm compared to other al- 
gorithms Ej for fixed point quantum search. It is also instructive to compare our results 
with a simple classical algorithm, where each iteration consists of picking an item randomly 
from the database and testing it by an oracle query, giving an error probability t q+l after q 
iterations. 

Real variables: No complex numbers appear in our algorithm, unless they are part of the 
operator U . This increases the ease as well as the possibilities for physical implementation 
of the algorithm. In particular, the algorithm can be realized using classical waves (with 
individual states represented as orthogonal wave modes), as is the case for the quantum 
search algorithm 0. 

Allowed values of q: A practical criterion for stopping the iterative algorithm would be 
that the error probability becomes smaller than some predetermined threshold e t h- Provided 
we have an upper bound, e < e up < 1, we can guarantee convergence by choosing e^ -1 " 1 < eth- 
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In our algorithm 2q + 1 can take all odd positive integer values, while in the Phase- 7r/3 search 
(or its alternate versions) 2q + 1 can only take a restricted set of values of the form 3 n (for 
positive integer n). Thus we can use q an — \\ [~ )°gg th — 111, while the Phase-7r/3 search 
requires q„ /3 = I(3r io g30°g Wi°g^ P )l _ X ) ( = 1,4,13,40,121,364,1093,...)- Thus we can 
save a sizable number of oracle queries — up to a factor of 3 — compared to the Phase- 7r/3 
search, especially when log 3 (logeth/loge up ) slightly exceeds an integer. Moreover, the simple 
classical algorithm needs q c \ — [ j"^" 1 — 11 oracle queries, which is always more than that for 
our algorithm (by about a factor of 2) but can be less than that for the Phase-7r/3 search (by 
up to a factor of 2/3). 

Deterministic vs. probabilistic: Although our complete algorithm is probabilistic, 
whenever we exit the iteration loop after finding ancilla-2 in the state |1), we obtain a deter- 
ministic result for the target state. The total probability of this deterministic result, over q 
iterations of our algorithm, is 



2fe-4 

k=2 



(8) 



Only the last measurement of the register gives a probabilistic result for the target state, with 
probability ie 2<? ~ 2 (l + e — 2e 3 ). The fact that a major fraction (asymptotically all) of the 
total 1 — e 2q+1 success probability gives a deterministic result can be a useful feature in some 
applications. Such a deterministic feature is not present in the Phase- 7r/3 search, but it exists 
in its alternate version 

Average number of oracle queries: The Phase-7r/3 search always requires the same 
number of queries as the number of iterations, because it has to go through all the iterations. 
On the other hand our algorithm can terminate at various times depending on the outcome 
of the ancilla-2 measurements. The average number of queries is always less than the number 
of iterations because our algorithm generally exits the iteration loop before completing it. We 
calculate the average number of oracle queries for our algorithm as in Eq.(8), 

^)+(^)<i-< 2 >x>"-*+^)^ 

l-e\ A + e\ A l-e 2 i- 2 \ . 1 - e 2 ?" 2 

(9) 



1-e 2 ) 2(1 -e) 



For the simple classical algorithm, the average number of oracle queries in case of 2q iterations 
is 

q cl = (1 - e) f] fee*" 1 + 2qe^ = ^! . (10) 
fe=l 

In general, 2q &n > q cll so our quantum algorithm can save at most a factor of 2 in the number 
of oracle queries compared to the simple classical algorithm. For q = 1, i.e. when the error 
probability is reduced to e 3 , g an = 1 is always better than q cl = 1 + e. For q > 1, our quantum 
algorithm provides an advantage only for e > e a , where e a satisfies 2e a + e 2<?_2 — 2e 2q — 1 and 
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Figure 4: Comparison of the average number of oracle queries for fixed point search algorithms, 
when q = 4. (a) is for our algorithm, (b) is for the Phase- 7r/3 search algorithm, and (c) is for 
the simple classical algorithm. 



is very close to but always less than 1/2. For e < e a , the simple classical algorithm requires 
less number of oracle queries than our quantum algorithm, but in this range q < 2 and hence 
the disadvantage is not significant. To illustrate this behavior, we have plotted the average 
number of oracle queries as a function of e in Fig. 4, for the three algorithms when q = 4. 

General ancilla-1 state: In the algorithm discussed so far, using ancilla-1 in the initial 
state |+), we effectively halved the fraction of target states. More generally, we can make 
the fraction of target states smaller by a factor r, by choosing ancilla-1 in the initial state 
y/l — r|0) + The corresponding joint unitary operator Uj becomes R®U : where R 

transforms |0) to \f\ — r|0) + \[r\\). It can be easily shown that after q iterations of the 
algorithm, the probability of getting a non-target state becomes e(l — 2r(l — e)) q . Also, the 
average number of oracle queries becomes (cf. Eq.(9)) 

1 - (1 - 2r(l - e)) 2 ?- 2 
1 - (1 - 2r(l - e)) 2 

= i + 1 -( 1 - 2 ;( 1 7^' 2 . (id 

4r(l-e) V ' 

Optimization of worst-case and average-case behavior: For a given initial error 
probability e, the optimal worst-case algorithm provides the fastest decrease of the error 
probability as a function of q, while the optimal average-case algorithm yields the smallest 
value of q. The above results give us the following solutions: 

(a) For r > 1/2, the error probability vanishes for a non-zero value of e, eo = 1 — V(2r)- So 
if e € (0, |] is known, a single oracle (q = 1) can find a target state with certainty. 

(b) For r > 1/2, there is a particular value of e, given by e T i = e (0, |], above which 



«a» - r(l-e) + (l-r(l-e)) 1 + 
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|1 — 2?'(1 — e)| < e. Thus, if we have a lower bound on e, choosing r to make it equal to e r i 
gives an algorithm that converges faster than e 2q+1 , i.e. the Phase-7r/3 search. It has been 
shown E] that for one iteration, the same behavior can be obtained using only one ancilla. 
In absence of a lower bound on e, our algorithm with r = 1/2 remains the best bet. 
(c) The optimal average case algorithm is easily obtained by taking the limit q — > oo. 

Clearly the smallest value for q is obtained for r — 1, which is the simple scheme without 
ancilla- 1 described at the beginning of section 2. It can give up to a factor of 4 improvement 
over the simple classical algorithm, when e > 1/4. The simple classical algorithm is better 
for e < 1/4, but in this range q < 4/3 and the difference is rather small. 

No. of ancilla states: Our algorithm uses only two ancilla qubits. The Phase-7r/3 
search needs a six-state ancilla to obtain the phase transformations from the standard binary 
oracle Of : \x)\a) — > \x)\a ®e f( x ))- Its alternative versions E] use one ancilla qubit per 
recursion step, and so require n ancilla qubits where 2q + 1 = 3™. Our algorithm is therefore 
quite economical in the number of ancilla states required. 

A related issue is whether or not the algorithm can be improved by POVM measurements. 
Since POVM measurements are equivalent to projective measurements in an enlarged Hilbert 
space (which can be obtained by adding extra ancilla qubits), and our algorithm has already 
achieved the optimal asymptotic performance with only two ancilla qubits, there is not much 
left for POVM measurements to improve up on. 

5 Outlook 

We have presented a new algorithm for fixed point quantum search using irreversible measure- 
ment operations, which is superior to the Phase-7r/3 search in several respects. Fixed point 
quantum search algorithms provide new techniques for driving a quantum state towards a 
specific target. Although they are no match for the optimal quantum search algorithm for 
large number of queries, i.e. 0(1/ f) compared to 0(1/ \fj) for small /, they are still better 
than the best classical search algorithm by a factor of 2. The potential applications of these 
algorithms are in problems with small but unknown initial error probability and limited num- 
ber of oracle queries, where they assure monotonic power-law convergence while the optimal 
quantum search algorithm is uncertain in improvement. This has been illustrated in case of 
correction of systematic quantum errors 0J|S]. 

The concept of ancilla-controlled quantum search presented here can be generalized to 
make the ancilla transformations iteration dependent, in a manner reminiscent of local adia- 
batic quantum algorithms, and that will be presented in a forthcoming paper 
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